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ABSTRACT 

We present the detection of four rotational emission lines of water vapor, from energy levels E u /k — 
101 - 454 K, in the gravitationally lensed z = 3.9 QSO host galaxy APM 08279+5255. While the 
lowest H2O lines are collisionally excited in clumps of warm, dense gas (density of hydrogen nuclei 
7iH = (3-1 ± 1.2) x 10 6 cm -3 , gas temperature T g ~ 105 ± 21 K), we find that the excitation of the 
higher lines is dominated by the intense local infrared radiation field. Since only collisionally excited 
emission contributes to gas cooling, we conclude that H2O is not a significant coolant of the warm 
molecular gas. Our excitation model requires the radiatively excited gas to be located in an extended 
region of high 100 fim opacity (tioo = 0.9 ± 0.2). Locally, such extended infrared-opaque regions are 
found only in the nuclei of ultraluminous infrared galaxies. We propose a model where the infrared- 
opaque circumnuclear cloud, which is penetrated by the X-ray radiation field of the QSO nucleus, 
contains clumps of massive star formation where the H2O emission originates. The radiation pressure 
from the intense local infrared radiation field exceeds the thermal gas pressure by about an order of 
magnitude, suggesting close to Eddington-limitcd star formation in these clumps. 
Subject headings: galaxies: ISM; galaxies: nuclei; quasars: individual (APM 08279+5255) 



1. INTRODUCTION 

Water is expected to be one of the most abundant 
molecules in molecular clouds in galaxies. In cold molec- 
ular clouds water is in the form of icy mantles on dust 
grains, with total H 2 a bundances up to 10 ~ 4 with re- 
spect to hydrogen nuclei (|Tielens et al.lll99lT ). thus con- 
taining up to 30% of the available oxygen atoms. In 
warm molecular clouds, such as in star forming galaxies 
or galaxies with a luminous active galactic nucleus, water 
can evaporate from the dust grains when the grain tem- 
perature becomes sufficiently high. H2O molecules can 
also be released into the gas phase by photodesorption 
in regions exposed to ultraviolet (UV) radiation and by 
desorption ind uced by cosmic rays or X -rays in more ob- 
scured regions (Hollcnb ach et al.ll2009D . or by sputtering 
of grains in shocks. In warm molecular gas, H2O can also 
be formed in the gas phase, through ion- neutral chem- 
istry in regions with a sufficiently high fractional ion- 
ization, or through neutral-neutral chemistry in regions 
sufficiently warm that the relevant activation energy bar- 
riers can be overcome. Gas-phase H2O may play an im- 
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portant role as a cooling agent of warm, dense molecula r 
clouds (|Neufeld fc Kaufman! Il993l : INeufeld et"aD[l995h . 
The large Einstein A- values of H2O rotational transi- 
tions lead to high critical densities (^ 10 s cm -3 ) so that 
collisional excitation will only be effective in very dense 
gas. 

Due to the wet Earth atmosphere, bulk gas-phase 
water can only be detected from space, or from dis- 
tant objects where the cosmological redshift moves the 
H2O lines into transparent atmospheric windows. De- 
spite the det ection of a 22 GHz water maser in a 
z = 2.6 QSO (jlmpel lizzeri et al. 2008), previous searches 
for non-maser rotational emission li nes of H 2 Q from 
high-z objects remained un successful ()Wagg et al.ll2006l : 
iRiechers et~aLll2006l l2009bh . until the recent detection of 
the 2o, 2 — li.i H2O emiss ion line from a gra vitationally 
lensed galaxy at z = 2.3 (jOmont et al.ll20lTI ). However, 
neither the molecular excitation mechanism nor the wa- 
ter abundance could be derived on the basis of the de- 
tection of a single optically thick line. 

The Herschel Space Observatory has recently enabled 
the first detections of H2O emission lines from two nearby 
galaxi es, with very different results between the two ob- 
jects (IVan der Werf et all 120101 : IGonzalez- Alfonso et al.l 
120101 : 1 Weifi et al.ll2010f ). Spectra of the lowest H 2 tran- 
sitions in the nearby starburst galaxy M82 revealed faint 
lines with complex spectral shapes, with one of the lines 
(the li.i ~0p n par a-HsO ground-state line) in absorption 
(jWeifi et al.ll2010i ). No H 2 Q emission from higher rota- 
tional levels was found (jPanuzzo et al.ll2010l ). In marked 
contrast, observations of the nearby Ultraluminous In- 
frared Galaxy (ULIRG) and QSO Mrk231 revealed a rich 
spectrum of H2O lines from upper level energies E n up 
to E u /k = 640 K (where k is the Boltzmann constant), 
whereas the lower lin es connecting to the grou nd state 
remained undetected (|Van der Werf et aLll2010f l. 
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Fig. 1.— H 2 emission from APM 08279+5255. The contoured 
pseudocolour map shows the distribution of the velocity-integrated 
continuum-subtracted H2O 2i,i — 2o,2 line flux, with contour lev- 
els of 3, 5, 7 and 9<r, where a = 0.5 Jy km s _1 beam -1 . The 
dashed white ellipse indicates the s ynthesized beam. The inset 
presents the NICMOS F110W image dlbata et aljl993 ) and shows 
that the gravitationally lensed images (the brightest of which are 
0'.'378 apart) are fully covered by the synthesized beam. 

In order to search for H 2 in a high-z galaxy in a 
systematic way, we undertook a search for four lines 
of H 2 in the ga s-rich z = 3.9 QSO host galax y 
APM 08279 + 5255 (|Irwin et alJll998tlLewis et al.lfl99l . 
This object is gravitationally lensed with a magnification 
factor i i — 4 according to the model of iRiechers et al.1 
(2009b). With CO lines detected from rotational levels 
up to J = 11 (jWeifi et al.l [20071 IRiechers et afll2009al ) 
and HC N, HNC and HCO+ lines from levels up to 
J = 6 (IWaggetal.1 [2005 IGarcia-Burillo et al.l [20061: 
IRiechers et al.l l2010f) . its molecular medium has been 
characterized very well. Analysis of the C O rotational 
ladder revealed unusually high excitation (Wei fi et al.1 
l2007h . However, a previous search for the ortho ground- 
state li,o — 1q,i H 2 l ine (upper level e nergy E u /k = 
61 K) was unsuccessful (I Wage et al.ll2006r i. We therefore 
targeted lines of higher excitation, with E u /k from 101 
to 454 K. 

2. OBSERVATIONS AND RESULTS 

We used the IRAM Plateau de Bure Interferometer 
(|Guilloteau et al.l 11992) with six antennas in December 
2010 and February 2011 to observe the four H2O lines 
listed in Table Q] (which also lists all relevant obser- 
vational parameters) in APM 08279+5255. Observing 
times (including overheads) varied from 2 to 5.2 hours 
per line. The WIDEX backend was used, providing 
3.6 GHz instantaneous bandwidth in dual polarization. 
Data reduction using the GILD AS package included the 
standard steps of data flagging, amplitude, phase and 
bandpass calibration, and conversion into datacubes. 
The flux density scale is accurate to within 15%. The 
datacubes w ere deconvolved using a CLEAN algorithm 
(IClarklll980n . 

All four lines were detected and an image of the flux 
distribution of the 2i i — 2o,2 line is shown in Fig. [TJ 
Fits to the continuum uw-data were used to calculate 
the continuum fluxes (reported in Table [JJ and source 
sizes. In the high angular resolution observations the 
emission was found to be slightly spatially resolved, with 
source sizes of approximately 0'.'5, in agreement with 
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Fig. 2.— Spectra of H2O emission lines from APM 08279+5255 
at a common frequency resolution of 40 MHz. The horizontal axis 
shows velocity relative to redshift z = 3.911. The red curves indi- 
cate Gaussian fits to each spectrum. 

high resolution o bservations of the 1 mm continuum 
(IKrfps et al.ll2007l) . and CO 1-0 (IRiechers et al.H2009bD . 
Integrated spectra of the four H 2 lines are shown in 
Fig. [5] and parameters of the detected lines are pre- 
sented in Table [TJ The most sensitive detection is that 
of the 2i 1 — 20,2 line, which displays a symmetric line 
profile with a full width at half maximum (FWHM) of 
480 ±58 km s -1 , excellently matching the FWHM value s 
of 445 - 520 km s" 1 of the CO lines (I Weifi et al.ll200l . 
Recently, a low spectral resolution 1 mm spectrum of 
APM 0812 79+5255 obtained wi th Z-Spec has been re- 
leased by Bradfor d et al.l (|201 If ) . In the Z-Spec spec- 
trum, the 32,i — 3i.2 line is the only H2O line detected 
at the 3cr level, but our flux meaurement for this line 
is a factor 2 lower. This discrepancy is similar to that 
between the CO line fluxes measured by Z-Spec and by 
the IRAM 30 m telescope, several of which have been 
independently confirmed by the IRAM Plateau de Bure 
Interferometer (Weifi et al. 2007). A similar discrepancy 
is found betwee n our flux measuremen t of the 42,2 — ^1.3 
line and that bv lBradford et al.l (|2011l) . although the lat- 
ter had a significance of only 2.5er. The continuum flux 
densities match very well between the two data sets. Fi- 
nally, we note that the serendipitous detection of another 
H2O line, the 22,0 — line, using the IRAM P lateau de 
Bure Interferometer, was recently reported by iLis et all 
(|20ll . 

3. FAR-IR PUMPED H 2 EMISSION 

The distribution of the detected line flux as a function 
of the energy of the upper level of the transition is shown 
in Fig. [3] The flux distribution shows considerable emis- 
sion out to the 4 2 ,2 — 4i,3 line, implying that for purely 
collisional excitation the kinetic temperature must be of 
the order of at least 400 K. Critical densities are ap- 
proximately 10 8 cm" 3 for the lowest lines (2n,2 — li.i 
and 2i.i — 20.2) and higher than 10 s cm -3 for the higher 
lines. Therefore, any purely collisionally excited model 
that can account for the observed 3 2 .i— 3i.2 and ^2,2—^1.3 
fluxes would produce much stronger emission in the lower 
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TABLE 1 

Parameters of H 2 line emission from APM 08279+5255 





2o,2 — ll.l 


2i,i — 2o,2 


32,1 — 3l,2 


42,2 — 4l,3 


li.o — lo.i 






987.926 


752.033 


1162.911 


1207.638 


556.936 


GHz 


b 

^obs 


201.166 


153.132 


236.797 


245.905 




GHz 


E u /k c 


101 


137 


305 


454 


61 


K 


Flux density calibrator 


3C273 


3C84 


3C273 


0923±392 






and flux density d 


8.0 


8.5 


7.3 


4.3 




Jy 


Synthesized beam c 


0'.'58 x 0'.'60, 33° 


2'.'98 x 1'.'87, 47° 


0'.'62 x 0'.'56, 48° 


0'.'61 x 0'.'44, 49° 






a { 


2.5 


1.7 


7.1 


4.3 




mjy 


Scont S 


16.5 ±0.8 


5.4 ± 0.3 


26.6 ± 1.3 


31.4 ±2.0 




mjy 


T h 

■"line 


9.1 ± 0.9 


4.3 ±0.6 


8.0 ±0.9 


7.5 ±2.1 


< 0.7 n 


Jy km s 


-^model 


8.9 


4.5 


8.4 


6.5 


0.25 


Jy km s 




72 ± 35 


47 ±25 


93 ±16 


50 ±43 




km s — 1 


FWHM k 


492 ± 84 


480 ± 58 


312 ± 37 


429 ± 96 




km s — 1 


L'y 

line 


(7.9 ±0.8) x 10 10 


(6.5 ± 1.0) x 10 10 


(5.1 ±0.5) x 10 10 


(4.5 ± 1.3) x 10 10 


< 1.9 x 10 10 


K km s" 


Mine 


(2.4 ±0.5) x 10 9 


(8.9 ± 1.3) x 10 8 


(2.6 ±0.3) x 10 9 


(2.6 ±0.7) x 10 9 


< 1.1 x 10 8 


Lq 



a Rest frequency of the line 

b Observing frequency 

c Upper level energy of the transition 

d Assumed flux density of flux calibrator at observing frequency 

c Full width at half maximum and position angle 

f R.m.s. noise in the integrated spectrum at Au = 40 MHz 

g Observed integrated continuum flux density 

h Observed integrated line flux 

1 Modeled integrated line flux 

J Center velocity of fitted Gaussian relative to z = 3.911 
k Full width at half maximum of fitted Gaussian 

I Using z = 3.911, H = 70 km s" 1 Mpc" 1 , n A = 0.73, Q t ot = 1 
m Not corrected fo r lensing 

II Upper limit from lWagg et al.l (2006) 
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Fig. 3. — Flux in H2O emission lines from APM 08279±5255 as 
a function of upper level energy E u . Black symbols and error bars 
indicate the observed line fluxes and the previously published upper 
limit on li,o — lo,l (indicated by a downward arrow). Red crosses 
indicate values produced by our model as described in the text. 
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tive upper limit exists (jWagg et al.i 120061 ) . 
fore rule out purely collisional excitation and consider in 
addition radiative excitation by an intense far-infrared 
(far-IR) radiation field from warm dust, with a temper- 
ature Td = 220 K, as d erived from the continuum spec- 
tral en ergy distribution ([Weifi et al.ll2007t iRiechers et al.l 
20093). 
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Fig. 4. — Partial H2O energy level diagram. The ortho and para- 
H2O rotational ladders are indicated, and the transitions detected 
here, as well as the upper limit on the li,o — lo,i ortho-H 2 ground 
state transition are indicated by red arrows. Purple arrows indicate 
the radiative pumping transitions at 58 and 75 fim that account 
for the efficient population of the 42,2 and 32,1 levels. 

We model th ese coupled processes usi n g a r adiative 
transfer code (jPoelman fc Spaansl 120051 120061 ) which 
computes the statistical equilibrium populations of all 
relevant H2O levels in the ground and first vibrationally 
excited state. Line radiation is transfcred in a non-local 
and three-dimensional manner through a multi-zone es- 
cape probability calculation. The model consists of dense 
clumps within a spherical region with a radius of 100 pc. 
Due to the high critical densities of the H2O lines, only 
the clumps contribute to the H2O emission, and we have 
run a model grid where we varied gas temperatures over 
the range T g = 50 — 200 K, and H2O column densi- 
ties over 10 16 - 10 18 cm . Clump densities were var- 
ied only over the range jih = 10 5 — 10 7 cm" 3 , since 
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higher densities (~ 10 8 cm' 3 ) would produce low rota- 
tional lines (21,1 — %o,2 and 2 ,2 — li,i) much stronger 
than 32,i — 3i,2, contrary to what is observed. The H2O 
ortho- to-para ratio was obtaine d from thermal equilib- 
rium (Poclman & Spaans 2005) and found to be very 
close to the statistical equilibrium value of 3. We fixed 
the gas-to-dust mass ratio at 100, and the local turbulent 
velocity dispersion at 4 km s _1 , with a total velocity dif- 
ference (gradient times length) across the model region 
of 140 kms" 1 . 

In our best fit model, the clumps have a density of 
hydrogen nuclei Tin = (3.1 ± 1.2) x 10 6 cm~ 3 , a gas tem- 
perature T g = 105 ± 21 K, and the average density of 
hydrogen nuclei over the 100 pc radius model region is 
(nii) = 4000 ± 1000 cm -3 , yielding a 100 /im continuum 
optical depth T100 ~ 0.9 ± 0.2. The line fluxes predicted 
by the best fit model are listed in Table [T] and shown in 
Fig. [3] The lowest lines (up to 2i,i — 2 ,2) are excited 
mostly by collisions. The clump density and tempera- 
ture derived in our model are determined by these lines, 
and the required density is somewhat higher than that 
derived from the CO lines (jWeifi et all 12007'). which is 
not surprising given the higher critical densities of the 
H2O lines. In contrast, the 42,2 and 32, 1 levels are popu- 
lated exclusively by the absorption of far-IR photons. As 
shown in Fig. 01 the pumping occurs at 58 and 75 rest 
frame wavelengths. The intensity of the radiatively ex- 
cited lines provides a measure of the intensity of the local 
far-IR radiation field and therefore of the optical depth 
in the spectral region of the pumping wavelengths, and 
this determines the 100 /im optical depth in our model. 
The relative strengths of the 4 2 ,2 — 4i, 3 and 32, 1 — 3i,2 
lines in principle provides information on the colour tem- 
perature of the far-IR radiation field, but for Td > 100 K 
the pumping wavelengths are on the Rayleigh- Jeans tail 
of the warm dust continuum, making the line ratio in- 
sensitive to the dust temperature. Since only collisionally 
excited lines contribute to the cooling (i.e., removal of ki- 
netic energy) of warm molecular gas, we find that cooling 
by the H2O lines is unimportant compared to the cool- 
in g by CO rotational line s, in contrast to the conclusion 
bv lBradford et all (120111) and to earlier theoretical sug - 
gcstions (INeufeld fc Kaufman|[l99l INeufeld et al.1119951) . 
We note that radiative excitation has also been suggested 
to drive the intensity of HCN lines in APM 8279+5255 
(|Garci'a-Burillo et all [20061: iWeifi et all [20071 ) . Further- 
more, the upper levels of the [Cn] 158 /im line and the 
[Si II] 35 /xm line are also affected by radiative pumping, 
as shown by rest frame UV abso rption line measurements 
by ISrianand k, Petitieanl (|2000D , who conclude that the 
absorbing clouds must be located closer to the extended 
(> 200 pc) source of far-IR radiation than to the UV 
source. 

In our best fit model, the 42,2 — 4i,3 line is only moder- 
ately optically thick (optical depth approximately 1.7), 
and is thus sensitive to the total column density of 
warm (T g > 100 K) water vapor in our model region 
(unlike the other lines, which have line center optical 
depths from 30 to 300). We find a total (beam-averaged) 
warm H2O column density of 2.6 x 10 17 cm~ 2 , i.e., a 
warm H2O/H2 adundance of approximately 6 x 10~ 7 , 
consistent with UV or X-r ay irradiated chemical models 
(jMeiierink & Spaansl [20051) . 



4. IMPLICATIONS FOR THE NUCLEAR REGION OF 
APM 08279+5255 

The H 2 line ratios observed in APM 08279+5255 
are very different from those observed in prominent 
star forming regions in the Milky Way. Water lines 
from UV irradiated gas (photon-dominated regions or 
PDRs ) in the Milky Way show therm al level popula- 
tions ([White et al.ll2010t lHabart et al.ll2010h . dominated 
by low-lying lines. These lines are much fainter than 
CO lines in the same frequency range. The best case in 
point is the prominent Orion Bar P DR where only a few 
low lying H 2 lines are detected ([Habart et alj l2010h . 
and the 2i,i — 2o, 2/11,0 — lo.i luminosity ratio is 0.6; in 
contrast, in APM 08279+5255, this ratio is > 8.1. Fur- 
thermore, in APM 08279+5255, the H 2 lines have lu- 
minosities comparable to CO lines in the same frequency 
range, while they are two orders of magnitude fainter 
in the Orion Bar (e.g., the H 2 2 M - 2 , 2 / CO(6-5) 
line luminosity ratio is 0.6 in APM 08279+5255, while 
it is 0.026 in the Orion Bar). The only object with 
properties comparable to APM 08279+5255 is the nearby 
ULIRG/QSO Mrk231 ([Van der Werf et alJlMol) . where 
high H2O lines were also foun d to be radiatively excited 
(|Gonzalez- Alfonso et al.|[2010D , and the two objects show 
very similar 32, 1 — 3i.2/2o,2 — I1.1 luminosity ratios. 

The fact that in the presence of an intense far-IR radi- 
ation field the most luminous H2O emission is radiatively 
excited follows directly from the large Einstein A- values 
of the H2O rotational lines, which result in high crit- 
ical densities as well as a strong coupling to the local 
far-IR radiation field. A key result of our analysis is 
therefore the presence of a sufficiently intense local far- 
IR radiation field, which requires the emitting clumps to 
be located in an obscured environment, with rioo ~ 0.9. 
The essential difference with the H2O excitation in the 
Orion Bar PDR (aside from less important differences in 
temperature and density) is the fact that the latter has 
a low 100 /im optical depth, while in APM 08279+5255 
Two ~ 0.9 over the entire region sampled. The min- 
imum size of the this region follows from the black- 
body limit, which is calculated by assigning the entire 
Lfir = 5 x 10 13 L (corrected for gravitational amplifi- 
cation) to a blackbody with T = 220 K. The resulting 
radius is R = 110 pc for a spherical model or a factor 
y/2 larger for a disk-l ike configuration, as preferred by 
iRiechers et al.1 (|2009bl ). The most recent lensing model, 
which accurately reproduces the positions and luminosi- 
ties of the gravitationally lensed images of the QSO nu- 
cleus, assigns the CO(l— 0) emission to a disk of radius 
of 550 pc, at an inclin ation of less than 30° from face-on 
([Riechers et al.ll2009b[ ). Because of the excellent match 
of the H2O and CO line widths, it is likely that the re- 
gion of H2O emission, with Td = 220 K and T100 ~ 0.9, 
has the same radius as the CO emission region. 

The CO rotational line ladder of APM 08279+5255 
shows even higher excitation than that of Mrk231 and 
therefore most likely reveals the presence of an XDR 
(|Weifi et al]l2H07l : IVan der Werf et al.ll27)Tol ). Since high 
H 2 a bundances can be generat ed in both PDRs and 
XDRs (IMeiierink fc Spaansl [2005). our data do not dis- 
tinguish directly between these models. However, while 
X-rays easily traverse and heat large column densities of 
gas, they are inefficient at heating the dust, and dust 
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temperatures higher th an ~ 50 K are hard to achieve 
over extended regions (jYanl 119971 iMeiierink fc Spaansl 
21 It ).->;•. The high dust temperature of ~ 220 K over a 
~ 550 pc radius region is more easily achieved in PDRs, 
generated by widespread star formation in clumps of 
dense gas throughout the circumnuclear gas disk. Re- 
gions of several 100 pc radius, with high 100 ^im optical 
depths, a re locally found only in the central regions of 
ULIRGs ([Solomon et all 119971 ). The fact that in this 
study we find such a region in a QSO, provides support 
for the scenario where U LIRGs form the birthplaces of 
QSOs (|Sanders et al.l[l988l) . 

If our suggestion that the extended warm dust con- 
tinuum is generated by circumnuclear star formation in 
APM 08279+5255 is correct, this star formation is taking 
place in the presence of the strong and p enetrative X-ray 
radia t ion field generated by the AGN ([G allaghe r et all 
[2001 IHasinger et all [2002t IChartas et all I2002D. As 
shown bv lHocuk fc Spaansl (I2010D . in this situation frag- 
mentation is inhibited and a top-heavy stellar Initial 
Mass Function is expected to result. It is possible 
that this effect accounts for the extraordinary far-IR 
luminosity Lfir = 5 x 10 13 L Q (corrected for gravi- 
tational amplification) of APM 08279+5255. We note 
also that in the clumps in our model, both turbulent 
pressure Pturb ~ P a tl^> (where p is the mass density 
of the gas clump and a v its turbulent velocity disper- 
sion) and radiation pressure P ra d ~ TioofT^/c (where 
a is the Stcfan-Boltzmann constant and c the speed 
of light) exceed the thermal pressure Pth ~ ri}i 2 kT g 
by a large factor. Inserting numbers from our best-fit 



model, we find that Pth ~ 3 x 10 -8 erg cm -3 , while 
Pturb ~ Prad ~ 3 x 10~ 7 erg cm -3 . Radiation pres- 
sure therefore plays an important role in the dynamics 
of the circumnuclear gas cloud. If the clumps are indeed 
forming stars, as we are suggesting, the star formation 
process is then close to Eddi ngton-limited, in agreem ent 
with the model developed by I Thompson et all (|2005l) for 
starburst regions surrounding a supermassive black hole. 

Our study demonstrates how radiatively exited H2O 
lines can be used to reveal the presence of extended 
infrared-opaque regions in galactic nuclei (even without 
spatially resolving these regions). Furthermore, we can 
derive local conditions in the infrared-opaque nuclear 
gas disk, which in the present case indicates close to 
Eddington-limited star formation. While for local galax- 
ies observations from space will remain necessary to ob- 
serve the H2O lines, the Atacama Large Millimeter Ar- 
ray will make this diagnostic readily available in galaxies 
with sufficient redshift, without the aid of gravitational 
lensing. 

This work is based on observations carried out with 
the IRAM Plateau de Bure Interferometer. IRAM is 
supported by INSU/CNRS (France), MPG (Germany) 
and IGN (Spain). DR acknowledges support from NASA 
through a Spitzer grant. We thank Melanie Krips for 
expert assistance with the IRAM data reduction, and 
Rodrigo Ibata for providing the NICMOS image used in 
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